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ABSTRACT 

The main theme of this paper concerns methods that may be classified as global 
(approximate) and local (exact). Some specific applications of these methods are 
found in: 

(1) Fracture and fatigue analysis of structures with 3-D surface flaws 

(2) Large-deformation, post-buckling analysis of large space trusses and space 
frames, and their control 

(3) Stresses around holes in composite laminates 
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A typical engineering problem is illustrated in figure 1, which shows a corner 
flaw at the intersection of a nozzle and a pressure vessel. The shape of the 
surface flaw may often be approximated mathematically as quarter-elliptical or 
quarter-circular. For the problem shown in figure 1, wherein the crack is located 
in the longitudinal plane of symmetry of the structure, only the so-called Mode I 
conditions prevail. In figure 1, the presence of a traction-free crack, in an 
otherwise unflawed solid, alters the stress-state only locally. From a viewpoint of 
fracture mechanics, however, the main quantities of interest are only the stress- 
intensity factors (strengths of asymptotic stress singularities) near the crack 
front. For analyzing fatigue crack growth and crack instability under thermal shock 
various flaw sizes and shapes need to be considered. The primary objective of 
analysis is to determine the variation of the Mode I stress-intensity factor along 
the border on the surface flaw. 




Figure 1. Corner surface-flaw at the pressure-vessel-nozzle intersection. 
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Figure 2 shows the schematic of a 12-bay space frame. The equations of dynamic 
motion of the frame, assuming large deformations and plasticity, may be written as: 


M (N+1 ) i + C (N+ 1 >X + 


C K 4X 


f -c + 


(N> r 




where M is the mass matrix, C the matrix of passive damping, t K the tangent stiff- 
ness matrix (which includes the effect of large deformation and plasticity), f^ the 
control-actuator force, Q„ the external load, the acceleration vector at 

time tj^ + 2 » v the ve l°city vector at AX the incremental displacement be- 

tween t N and t N+ j, and ^ 'R the internal-force vector at t N . In order to implement 
the control algorithms in an efficient manner, the order of the above system of 
equations must be as small as possible (i.e., each frame member must be modeled by 
no more than one finite element). Further, the control must be implemented for 
pulse-type loading of high intensity, such that the above system of equations must 
be integrated directly rather than using a modal-decomposition. Also requirements 
of on-line control may necessitate that C K, C, and M be known explicitly (in closed 
form) for arbitrary values of deformation, without the need for introducing approxi- 
mate shape functions for deformation of each element and without the need for any 
numerical integrations over each element. In figure 2, the object of inquiry is 
what effect does local (member) instability have on global (system) stability? How 
can we control the dynamic deformations locally to improve global behavior? Each 
member may be treated as a truss member, or a 3-D beam-type member, depending on 
joint design. How can local effects be accounted for simply and efficiently, so 
that algorithms for control of dynamic motion may be implemented, on line, using 
on-board computers in a large space structure? 
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An appraisal of the computational mechanics methods is given in figure 3. 
These methods include the f inite- element , boundary - element , and edge function 
methods (fig. 3). 


Finite Elements: 

• Trial and test functions are both approximate 

• Trial and test functions are, in general, alike - Galerkin 

APPROACH 

• In some instances it is best to have test functions different 

FROM TRIAL FUNCTIONS - PeTROV~GaLERKI N APPROACH 

• THE SOLUTION IS BOTH GLOBALLY AND LOCALLY APPROXIMATE 

• Versatile or arbitrary geometry, boundary conditions, suited 

FOR GLOBALLY approximate nonlinear solutions 

Boundary Elements: 

• Test functions are globally exact for the given linear problem, 

OR AT LEAST FOR THE HIGHEST'ORDER DIFFERENTIAL OPERATOR OF 
THE PROBLEM 

• Trial functions are approximate (at boundary only for linear 

PROBLEMS, AND IN INTERIOR ALSO FOR NONLINEAR PROBLEMS) 

• The SOLUTION IS BOTH locally and globally approximate 

• Not as versatile as the finite-element method, but excellent 

FOR SOME SPECIFIC PROBLEMS 

Edge Function Method: 

• Trial functions are globally exact 

• Test functions assumed only at boundary 

• Limited to linear problems posed by classical differential 

EQUATIONS 


Figure 3. Appraisal of computational mechanics methods. 
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In the most commonly used Galerkin finite-element approach in computational 
solid mechanics, the trial and test function spaces are identical and consist of 
simple piecewise continuous algebraic polynomials over each finite element, such 
that these functions and their appropriate-order derivatives (as dictated by the 
problem on hand) are continuous at the interelement boundaries. For problems of 
fourth or higher order, such as those of plates and shells, the development of 
finite elements has long been, and continues to be, plagued by the need for C 1 
(or higher order) continuity at the intereleraent boundaries. However, the success 
of the finite-element method in structural mechanics is unparalleled and is mainly 
due to the intuitive and ’geometric 1 interpretation of the method. The method is 
versatile in its ability to deal with complicated structural assemblies, such as of 
beams, plates, and shells, of the type used in aerospace applications. The solu- 
tions obtained through the finite-element method may be classified, in general, as 
being both globally as well as locally approximate. 

On the other hand, in linear and nonlinear solid mechanics, it is often pos- 
sible to derive certain integral representations for displacements. A key ingre- 
dient which makes such derivations possible is the singular solution, in an infinite 
space, of the corresponding differential equation (in certain linear problems) or of 
the highest-order differential operator (in the nonlinear case, or even in the 
linear case when the full linear equation cannot be conveniently solved), for a 
’unit* load applied at a generic point in the infinite space. When the problem is 
linear and the singular solution can be established for the complete linear dif- 
ferential equation of the problem, the aforementioned integral representations for 
displacements involve only boundary integrals of unknown trial functions and their 
appropriate derivatives. Such an integral representation, when discretized, leads 
to the so-called boundary-element method. Such pure boundary-element methods are 
possible in linear, isotropic, elastostatics , and in problems of static bending of 
linear elastic isotropic plates. On the other hand, as in such cases as (i) linear 
problems wherein the singular solutions cannot be established for the entire dif- 
ferential equations, (ii) anisotropic materials, and (iii) problems of large 
deformation and material inelasticity, the integral representations (if any) for 
displacements would involve not only boundary integrals but also interior-domain 
integrals of the trial functions and/or their derivatives. A discretization of such 
integral equations would lead not only to a simple boundary-element method but also 
to a sort of hybrid boundary/interior element method. 

When asymptotic solutions to the governing differential equations of the prob- 
lem are used as assumed trial functions, the interior residual error is zero, and 
only the boundary conditions need to be satisfied in a weighted residual method. 

Such an approach is called the edge-function method, but is limited mostly to linear 
problems. For further details, see references 1 through 4. 
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The surface-flow problems for current methods are noted in figure 4. 
for the proposed method are also shown. 


All Present Methods : Singular stress-state near the flaw 

border is modeled by locally approximate methods 


(I) Finite-Element Methods (singular elements) 
fttluri & Kathiresan (refs. 5-10) 

(Hybrid crack elements) 3-6,000 d.o.f. 
Tracy, Barsoum, Newman & Raju (refs. 11-13) 
(Distorted isoparametric elements and 
singular shape fn.) 5-10,000 d.o.f. 

These are very expensive, but accommodate 
arbitrary geometries of structure and flaw. 

(II) Boundary-Element Methods (for locally approxi- 
mate stress analysis and K-estimation from 
stress extrapolation) 

Cruse (ref. 14), Heliot et al . (ref. 15). 

Not suitable for 'thin' shells with flaws. 
Still very expensive. 

(Ill) Line-Spring Method 

Limited" to simple" geometries of structure and 
flaw. 


Prop o s ed Method : 

• It is a GLOBALLY APPROXIMATE, but LOCALLY EXACT 

METHOD 

• Singular stress-state near the flaw is NOT MODELED 

NUMERICALLY 

• It is about 30 times cheaper than the singular 

finite-element method 

• Details (Atluri & Nishioka (refs. 16-21) - several 

papers with varied examples) 


Figure 4. Surface- flaw problems. 


Problems 
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LLU 



Locally exact solution 
for pressurized crack 


Remarks: 1. Solution D: A rather complicated analytical 

solution (Atluri & Vijayakumar, Journal of 
Applied Mechanics , 1981) (ref. 2D 

2. Local solution due to crack-face traction alone 
is (i.e., the Solution C) the source of singu- 
larity. The stresses due to this local solution 
decay very rapidly. Only one or two i terations 
are sufficient to obtain K-solutions with 1%~ 
accuracy . 

About 30 times cheaper than the usual finite- 
element method for a typical problem such as 
flawed BWR nozzle. 


Figure 5. Global (approximate) and local (exact) analyses of embedded flaws. 
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Some comments concerning the solution of surface flaws in finite bodies using 
the present procedure are in order (fig. 4). Since the analytical solution of an 
elliptical flaw embedded in an infinite solid is used as solution D, it is necessary 
to define the residual stresses over the entire crack-plane including the fictitious 
portion of the crack which lies outside of the finite body containing only a surface 
flaw (fig. 5). Moreover it is well known that the accuracy of the 'least-squares' 
type function-interpolation inside the interpolated region can be increased with the 
number of polynomial terms; however, the interpolating curve may change drastically 
outside of the region of interpolation. The optimum variation of pressure on the 
crack surface extended into the fictitious region should be as shown in figure 6. 

For in-depth discussions of a variety of surface problems and their solutions, see 
refs. 19-21, 23-28. 



Figure 6. Postulated residual stress distributions on fictitious portions of an 
elliptical flaw, in the case of semi-elliptical or quarter-elliptical 
surface flaws. 
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For both 3-dimensional truss and frame members, explicit (locally exact) 
tangent stiffness matrices have been derived (fig. 7). Some effects of local 
(member) buckling on global (structural) behavior are illustrated in figures 8 
and 9. 


Truss Member: 

s ”Tach member undergoes large displacement and large 
rigid rotation 

• Member material is nonlinear 

• Each member may buckle and become curved (what 

effect does it have on global stability?) 


Frame Member: 

"Concepts of 3-dimensional semi-tangential rotations 
employed 

• Each member undergoes arbitrarily large rigid rota- 

tions and rigid displacements 

• Bending-stretching coupling incorporated in each 

member 

• Plastic-hinge method used to account for plasticity 

in each member 

• Member forces: axial, shear, and bending- twisting 

moments 


Kondoh and Atluri (refs. 29-30) and Kondoh et al . (refs. 31-32) 


Figure 7. Space trusses and space frames. 


Ill 




L = 66 04 (cm) 

(a) Thompson's Strut 
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(a) Load System: (i) P^: Vertical (b) Stability boundary under loads 

point loads as all nodes; and P 2 

(ii) P^: Vertical point loads at 
nodes in quadrants xj, x 2 > 0 


Figure 9. Study of the effect of member buckling on global (system) stability. 
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Examples of the efficiency of the global/local approach in analyzing frames are 
illustrated in figures 10 and 11. In figure 10, the classical problem of a two-bar 
frame is schematically illustrated. In the present approach, the tangent stiffness 
matrix of each member (respresented by a single finite element) is derived from 
exact solutions of governing differential equations which account for the bending- 
stretching coupling. Thus, no "shape functions 11 are assumed in each element, and no 
numerical integrations are performed in forming the tangent stiffness matrix. The 
present numerical integrations are performed in forming the tangent stiffness 
matrix. The present numerical results are shown to agree excellently with those of 
Wood and Zienkiewicz (ref. 34), as well as the experimental results of Williams 
(ref. 35). However, Wood and Zienkiewicz use five finite elements to model each 
member of the frame. 


In figure 11, the problem of plastic collapse of a frame is illustrated. Here 
again, the tangent stiffness matrix of each member (represented by a single finite 
element) is derived in closed form, accounting for large deformations, large rota- 
tions, and plasticity. A plastic-hinge method is used, and the progressive develop- 
ment of plastic hinges, at various load levels, is indicated in figure 11. 


^ai0243(,n) 



P(ib) 



Figure 10. Variation of load-point displacement and support reaction with applied 
load in a two-bar frame. Tangent stiffness of each 3-D beam member 
undergoing large deformation, large rotation, and plasticity is exact. 
Locally exact solution (ref. 30). 
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Figure 12 shows a problem of current interest in the analysis of stiffened 
composite plates. Issues involve the following: (1) stress concentrations near the 

hole in a composite laminate, (2) local buckling of stiffeners, (3) effect of geo- 
metric imperfections, (4) effect of discontinuities, and (5) three-dimensional 
effects and delaminations near the hole. An efficient globally approximate and 
locally exact approach could possibly include: (1) use of locally exact, laminated 

hole elements with embedded three-dimensional stress state (refs. 36 and 37), 

(2) use of locally exact stiffener elements as described earlier (ref. 32), 

(3) techniques for proper interacting of various elements, and (4) hole elements 
that can be improved by incorporating possible free-edge singularities in °^± m 


Focus Problem 



(GLOBALLY & LOCALLY APPROXIMATE) FINITE ELEMENTS] 

Figure 12. A stiffened laminated-composite panel with a hole. 
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Another example of the advantages of using a global/local approach is illus- 
trated here in the problem of analysis of stresses near a hole in a laminated 
composite [two cases of (-45/+45) g and (90/0 ) g laminates are discussed]. Figure 13a 
shows a typical finite-element model with "special-hole elements" in which a 3-D 
asymptotic hole solution is embedded. Figures 13b and 13c illustrate the excellent 
accuracy obtained from the present approach, in comparison with a fully 3-D finite- 
element solution of Rybicki and Hopper (ref. 38). The present solution is, however, 
an order of magnitude less expensive. Further details are given in references 36 
and 37. 



O' 


(a) (b) (c) 

Typical FlM model of a laminate Stress around a hole Stress around a hole 

with hole. 3-D asymptotic "hole- in C-45/+45) s laminate in (90/0) s laminate 

solutions" embedded in elements 
NEAR THE HOLD (Ref. 36) 


Figure 13. Analysis of stress state near a hole in laminated composites. 
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The following conclusions and recommendations are given. 

* Hybrid analytical/numerical methodologies should be explored 

* Simplified analysis procedures for elasto-plastic should be considered 
(Dynamic response calculations should be studied (some benchmark problems 
essential) ) 

* Constitutive models badly need improvement 

* Methods of coupling of problem-specific methodologies for use in general 
purpose programs should be explored 

* Trends to treat structural mechanics problems as continuum mechanics 
problems should be critically reviewed; the knowledge base in structural 
mechanics should be fruitfully utilized 

* Attempts to bridge the gap between micromechanics and macromechanics of 
heterogeneous (composite) media through computational mechanics should be 
explored 

* Computational stochastic structural analysis methods should be developed 

* Algorithms for new computing systems (MIMD) should be explored 

* Expert systems, . . . . (?) 

NASA’s role should be to provide: 

* Increased research support to universities 

* Predoctoral NASA fellowships (up to 20K per year, tax-free) that could be 
awarded to attract the best students 

* Long-range funding to properly plan and sustain high-quality research 
efforts 

* Increased access to supercomputers 

* Frequent visits to NASA facilities by graduate students to participate in 
laboratory testing. University facilities in this area are scarce; 
students in computational mechanics should get some first-hand experience 
in experimental mechanics 
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